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^ We present preliminary results of 7 years of Arecibo timing of the pulsar-white dwarf 

• , binary PSR J1738+0333. Wc can measure the proper motion, parallax with excellent 

' precision and have detected the orbital decay. Furthermore, the companion has been 

detected at optical wavelengths and a mass ratio of 8.1 ± 0.3 has been measured from 
the orbital variation of its Doppler shift. Once the companion mass is determined from 
< I , the optical measurements, this system will provide strong limits for the radiation of 

' dipolar gravitational waves. Assuming that general relativity holds, the fast-improving 

^ ' measurement of the orbital decay, combined with the measurement of the mass ratio, 

will provide an independent and precise measurement of the component masses. 
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\^ • PSR J1738-f0333 is a 5.85 ms pulsar in a low eccentricity binary system with or- 

, bital period of P(, = 8.5 iP and a low mass helium white dwarf (WD) companion. 

The pulsar was discovered in 2001 in a 20-cm intermediate Galactic latitude sur- 
vejOl carried out with the multi-beam system of the Parkes Telescope. Given its 
low Northern declination, it can be observed with much improved sensitivity using 
the the 305-m Arecibo Telescope. For that reason it has been timed with Arecibo 
^ . since 2003 using the L-wide receiver (with frequency coverage from 1100 to 1700 

MHz). We used the Wide-band Arecibo Pulsar Processor (WAPP) correlators and 
the Arecibo Signal Processor (a coherent dedispersion machine) as back-ends. The 
WAPPs have bandwidths of 100 MHz centered at 1170, 1410 and 1510 MHz (and, 
occasionally, a fourth WAPP centered at 1310 MHz as well). 

The WAPP correlator data were Fourier transformed to obtain one 256-point 
power spectrum every 64 /is. These are then folded and dedispersed at the pulsar's 
spin period and DM using "sigfoldpow" , a program designed by Ingrid Stairs specifi- 
cally for producing highly accurate average pulse profiles from search data. The ASP 
data have not yet been fully analyzed, that is one of the reasons why these results 
are still preliminary. The observed pulse profiles are then cross-correlated with a 
low-noise template, the best fit yields a pulse time of arrival at the telescope (TOA). 



1. Introduction 
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The TOAs were then analysed using the TEMPCB timing software, which makes a 
least-squared fit of the observed barycentric TOAs to a timing model, which includes 
spin, astrometric and orbital parameters. 

2. Results 

We have obtained a TOA residual rms of 1.4 /is, with TOAs taken every 268 
s for each of the 12 (or 16) simultaneous 12.5-MHz sub-bands. This results in 
more than 16000 TOAs. Of these, we use 9000 with a stated accuracy better 
than 3 /its. This large number of precise TOAs provide extraordinarily precise tim- 
ing parameters. These include the proper motion {fia = 7.041 ± 0.007 mas yr~^ 
and Us — 5.135 ± 0.019 mas yr~^) and parallax (tt = 0.75 ± 0.05 mas), which cor- 
responds to a distance oi d — 1.33 ± 0.09 kpc. The orbital period is measured 
so precisely {Pb = 30653.9194513 ± 0.0000005 s) that we can aheady detect the 
orbital decay, albeit to only 2-a precision: Pt — (—2.1 ± 1.0) x 10~^*ss~^. 
This becomes more significant after subtracting the extrinsic orbital period vari- 
ation caused by the proper motion of the system^ and the Galactic accelera- 
tion: Pb = (0.73 ± 0.08) X lO^^^ss^^. The intrinsic orbital period derivative 
is A'"* = (-2.9 ±1.0) X 10-14 ss-i. 

3. Optical Observations 

The companion white dwarf has been detected, and it has deep spectral lines typ- 
ical of a 0.2 M© thick-envelope He white dwarf. These were used to determine 
a radial-velocity curve using Gemini South on Cerro Pachon. From this we can 
derive the mass ratio of the system, R — 8.1 ± 0.3 (Van Kerkwijk, 2007, pers. 
comm.). At the moment we don't have a good estimate of the companion mass. 
However, if it follows the relation postulated by Tauris & Savonije in 1999,^^ then 
the companion mass rric — 0.18 ± 0.02 M©, which is qualitatively consistent with 
the observed spectral properties. Given the mass ratio this results in a pulsar mass 
Trip ~ 1.46 ± 0.16 A/q and an orbital inclination i ~ 33°. Apart from resulting in 
a reasonable value for irip, the low inclination is consistent with the non-detection 
of the Shapiro delay. Furthermore, when we introduce Shapiro delay parameters 
rric and sini in our solution, the apparent eccentricity diminishes to an insignifi- 
cant level of 0.00000039 ± 0.00000013. This implies that the difference between the 
semi-major and semi-minor axes is only 8 ± 6 ^m! 

4. Limit on dipolar gravitational wave emission 

If these values for nic and rrip hold, the orbital decay due to the emission of 

• GR 

quadrupolar gravitational waves predicted by general relativity is = — 2.7 x 

10"^'* s s""'^. This is consistent with the intrinsic orbital decay Pb 
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Fig. 1. Limits on G and Krt derived from the joint analysis of the timing data for PSR J0437— 4715 
plus PSR J1012+5307 (thin contours) or PSR J1738+0333 (thick contours). The inner contours 
include 68.3% and the outer contours include 95.4% of the simulated probability in each case. 



Subtracting Pb from Pb we obtain an "excess" of Pb — (—0.2 ± 1.0) x 
lO^^'^ss"^. Following the logic of Lazaridis et al. (these proceedings), we can 
combine this "excess" with that from a pulsar with a longer orbital period, like 
PSR J0437— 4715,'^ and derive joint limits for the variation of the gravitational 
constant G and the dipolar gravitational wave emission constant kd. In fig- 
ure [U we can see that, provided that the companion mass estimate is realistic, 
PSR J1738+0333 already provides excellent limits on dipolar gravitational wave 
emission: kd = l-^t^ij x 10"^ (1-cr) and kd = '^Atf^^ x 10""* (2-cr). This hap- 
pens despite it having a baseline that is only half that of PSR J1012+5307: this is 
more than compensated by its shorter orbital period and better timing precision. 
This also means that the measurement of Pb will keep improving much faster 
for PSR J1738-I-0333 in the future. One of the consequences is that, if we assume 

• int • GR 

general relativity to hold (i.e., Pb — Pb ), we will obtain precise values for rric, 
TUp and i from the measurements of Pb and R. These will help to calibrate the 
spectral models of low-mass WDs. 
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